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Interpretation of extensive air showers (EAS) experiments results is strongly based on air shower simulations. 
The latter being based on hadronic interaction models, any new model can help for the understanding of the nature 
of cosmic rays. The EPOS model reproducing all major results of existing accelerator data (including detailed data 
of RHIC experiments) has been introduced in air shower simulation programs CORSIKA and CONEX few years 
ago. The new EPOS 1.99 has recently been updated taking into account the problem seen in EAS development 
using EPOS 1.61. We will show in details the relationship between some EPOS hadronic properties and EAS 
development, as well as the consequences on the model and finally on cosmic ray analysis. 



1. INTRODUCTION 

Air shower simulations are a very powerful 
tool to interpret ground based cosmic ray experi- 
ments. However, most simulations are still based 
on hadronic interaction models being more than 
10 years old. Much has been learned since, in par- 
ticular due to new data available from the SPS 
and RHIC accelerators. 

In this paper, we discuss air shower simula- 
tions based on EPOS, the latter one being a 
hadronic interaction model, which does very well 
compared to RHIC data [T], and also other par- 
ticle physic experiments (especially SPS experi- 
ments at CERN). But used in air shower simula- 
tion program like CORSIKA [2] or CONEX [3], 
some results where in contradiction with KAS- 
CADE data [3], while it was better for other ex- 
periments [S]. 

Due to the constrains of particle physics, air 
shower simulations using EPOS present a larger 
number of muons at ground [6]. On the other 
hand, we will explain in this paper, how the 
contrains given by cosmic ray experiments can 
compensate the lack of accelerator data in some 
given kinematic regions (very forward) to improve 
hadronic interaction models and in particular the 
new EPOS 1.99. 



2. EPOS Model 

One may consider the simple parton model to 
be the basis of high energy hadron-hadron inter- 
action models, which can be seen as an exchange 
of a "parton ladder" between the two hadrons. 
In EPOS, the term "parton ladder" is actually 
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Figure 1. Elementary parton-parton scattering: 
the hard scattering in the middle is preceded by 
parton emissions attached to remnants. The rem- 
nants arc an important source of particle produc- 
tion even at RHIC energies. 

meant to contain two parts [7]: the hard one, as 
discussed above, and a soft one, which is a purely 
phenomcnological object, parameterized in Regge 
pole fashion. 
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In additions to the parton ladder, there is 
another source of particle production: the two 
ofF-shell remnants, see fig. [T] We showed in ref. 
[5] that this "three object picture" can solve the 
"multi-strange baryon problem" of conventional 
high energy models, see ref. [9]. 

Hence EPOS is a consistent quantum mechan- 
ical multiple scattering approach based on par- 
tons and strings [7], where cross sections and the 
particle production are calculated consistently, 
taking into account energy conservation in both 
cases (unlike other models where energy conser- 
vation is not considered for cross section calcu- 
lations [lO]). Nuclear effects related to Cronin 
transverse momentum broadening, parton satu- 
ration, and screening have been introduced into 
EPOS [11]. Furthermore, high density effects 
leading to collective behavior in heavy ion col- 
lisions are also taken into account [12] . 
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Figure 2. Total cross section of proton-carbon in- 
teractions. EPOS 1.99, QGSJETII, QGSJETOl 
and SIBYLL 2.1 hadronic interaction models 
(lines) are compared to data [13] (points) 

Energy momentum sharing and remnant treat- 
ment are the key points of the model concerning 
air shower simulations because they directly in- 
fluence the multiplicity and the inelasticity of the 
model. At very high energies or high densities, 
the so-called non-linear effects described in [TT] 
are particularly important for the extrapolation 
for EAS and it's one of the parts which has been 
changed in EPOS 1.99. 



2.1. Cross section and inelasticity 

We learned from KASCADE data 4 , that the 
energy carried by hadrons in EPOS 1.61 simula- 
tions is too low. It means than the showers are 
too old when they reach ground and it was due to 
a problem in the calculation of the nuclear cross 
section and to a too large remnant break-up at 
high energy (leading to a high inelasticity). 

To improve the predictive power of the model, 
the effective treatment of non-linear effects de- 
scribe in [11 has been made consistent to describe 
both proton-proton, hadron-nuclcus and nucleus- 
nucleus data with a unique saturation scale which 
can be fixed thanks to proton-proton cross section 
and Cronin effect in dAu collisions at RHIC. De- 
tails will be published in a dedicated article. 
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Figure 3. Inelastic cross section of proton-air in- 
teractions. EPOS 1.99, QGSJETII, EPOS 1.61 
and SIBYLL 2.1 hadronic interaction models 
(lines) are compared to data of air shower exper- 
iment (points). 

The EPOS 1.99 (full line) proton-carbon to- 
tal cross section is shown Fig [2| It is now 
in very good agreement with the data 13J and 
with the other hadronic interaction model used 
for air shower physics QGSJETOl [H] (dash- 
dotted line), QGSJETII [15] (dashed line) and 
SIBYLL [11] (dotted line). In fig d the ex- 
trapolation to proton-air data up to the high- 
est energies is shown in comparison with mea- 
surement from cosmic ray experiments. The er- 
ror bar for EPOS 1.99 represents the uncertainty 
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due to the definition of the inelastic cross sec- 
tion as measured by cosmic ray experiment^ In 
comparison with EPOS 1.61 (dash-dotted Hne), 
the EPOS 1.99 cross section has been notably re- 
duced. 

2.2. Particle production 

Thanks to a Monte Carlo, first the collision 
configuration is determined: i.e. the number of 
each type of Pomerons exchanged between the 
projectile and target is fixed and the initial en- 
ergy is shared between the Pomerons and the two 
remnants. Then particle production is accounted 
from two kinds of sources, remnant decay and cut 
Pomeron. A Pomeron may be regarded as a two- 
layer (soft) parton ladder attached to projectile 
and target remnants through its two legs. Each 
leg is a color singlet, of type qq , qqq or qqq from 
the sea, and then each cut Pomeron is regarded 
as two strings, cf. Fig. |4^. 
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Figure 5. Model comparison: longitudinal mo- 
mentum distributions of pion carbon collisions 
at 100 GeV from EPOS with (full) or without 
(dash-dotted) sting-end diquarks and QGSJETII 
(dashed) compared to data [TT] . 
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Figure 4. a) Each cut Pomeron is regarded as 
two strings b). c) The most simple and frequent 
collision configuration has two remnants and only 
one cut Pomeron represented by two q— q strings, 
d) One of the q string ends can be replaced by a 
qq string end. e) With the same probability, one 
of the q string ends can be replaced by a qq string 
end. 



'^The diflference between the top and the bottom of the 
error bars is the part of the cross-section where secondary 
particles are produced without changing the projectile 
(target dilfraction). Cross section of other models includes 
this target dilfraction (top of error bars). 



It is a natural idea to take quarks and an- 
tiquarks from the sea as string ends for soft 
Pomeron in EPOS, because an arbitary number 
of Pomerons may be involved. In addition to this 
soft Pomerons, hard and semihard Pomerons are 
treated differently. 

Thus, besides the three valence quarks, each 
remnant has additionally quarks and antiquarks 
to compensate the flavours of the string ends, 
as shown in fig. According to its number of 
quarks and antiquarks, to the phase space, and to 
an excitation probability, a remnant decays into 
mesons, (anti)baryons f8\. Furthermore, this pro- 
cess leads to a baryon stopping phenomenon in 
which the baryon number can be transfered from 
the remnant to the string ends (for instance inUji, 
depending on the process, the 3q-|-3q can be seen 
as 3 mesons or a baryon-antibaryon pair). 

In case of meson projectile, this kind of diquark 
pair production at the string ends leads to an in- 
crease of the (anti)baryon production in the for- 
ward production in agreement with low energy 
pion-nucleus data [lT| as shown fig. [5l As a con- 
sequence it is part of the larger number of muons 
in EAS simulations with EPOS. 

Compared to EPOS 1.61, EPOS 1.99 has a re- 
duced excitation probability at high energy, in- 
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creasing the number of protons in the forward 
direction and reducing the inelasticity. 

3. AIR SHOWERS 

In the foUowing, we discuss air shower simu- 
lations, based on the shower programs CONEX, 
using the old EPOS 1.61 (dash-dotted line), the 
new EPOS 1.99 (full hne) and QGSJETII(dashed 
line) (as a reference) as high energy hadronic in- 
teraction model in KASCADE experiment energy 
range. 
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Figure 6. Maximum hadron energy as a function 
of the primary energy for proton induced showers 
using EPOS 1.99 (fuU line), EPOS 1.61 (dash- 
dotted line) and QGSJETII (dashed hne) as high 
energy hadronic interaction model. 

The effect of the reduced cross section and in- 
elasticity is clearly visible on the maximum en- 
ergy of hadrons at ground as shown fig. [H The 
shower being younger at ground with EPOS 1.99, 
the maximum energy is up to 60% higher than in 
the previous release 1.61. The results are now 
close to QGSJETII results but with a diff'erent 
slope due to a different elongation rate. As a 
consequence, EPOS 1.99 does not have the prob- 
lems pointed out in [3] anymore and should be 
compatible with KASCADE data. 

4. Summary 

EPOS is a new interaction model constructed 
on a solid theoretical basis. It has been tested 
very carefully against all existing hadronic data. 



also those usually not considered important for 
cosmic rays. In EAS simulations, EPOS provides 
more muons than other models, which was found 
to be linked to an increased diquark production in 
both string ends and string fragmentation. The 
new EPOS 1.99 has a reduced cross section and 
inelasticity compared to the previous EPOS 1.61 
which leads to deeper shower development. This 
would solve the problem with KASCADE data. 
But since the number of muons and the elonga- 
tion rate are different than in the other models, 
the resulting analysis will be significantly differ- 
ent. 
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